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Abstract 

The distribution of radiation is investigated for the modeless laser 
having a multilobe mirror with the lobes (planes) inclined by small 
angles to optical axis. It is shown that change of the direction resulting 
from many passages of a ray through the optical system including 
a multilobe mirror may be described as Brownian walk of a point 
along the plane or equivalently as a solution of the two-dimensional 
diffusion equation. Boundary conditions for the diffusion equation 
may be approximately formulated as null conditions at some angle 
which, if being reached during the walk, guarantees that the ray escape 
from the optical system. In the framework of this approximation an 
explicit formula for the distribution of the outgoing ray in different 
angles is derived. After many passages through the optical system the 
angular distribution tends to some universal function. In the case of 
the round mirror it may be presented by the Bessel function of order 
zero. 
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1 Introduction 



It is known |Jl| that the angular distribution of the radiation of the laser with 
two plane-parallel mirrors is nonregular (inhomogeneous) because of nonideal 
mirrors and inhomogeneity of the active medium. This defect is absent in so 
called modeless lasersO having no longitudinal modes structure. Another 
configuration of modeless laser was proposed in the work 0]. In such a laser 
one of the plane mirrors is replaced by an multilobe mirror having the lobes 
(planes) slightly inchned to the optical axis (Fig.|T]a). It was experimentally 
shown 1^ that this replacement makes the angular distribution of the laser 
radiation smooth (homogeneous). There is no longitudinal modes in this laser 
too. Light is mixed and scattered in the cavity of a multilobe mirror laser. 

The form of the angular distribution in the multilobe laser was analyzed 
in the work by consideration of geometry of the rays crossing the optical 
system with the multilobe mirror many times. Such an analysis proves to be 
complicated, may be carried out only for some simplest configurations of the 
multilobe mirror and gives no explicit formula for the angular distribution. In 
the present paper we shall analyze the case of an arbitrary multilobe mirror 
and calculate an approximate angular distribution in an explicit analytical 
form. The idea of the method is following. 

In each double passage through the optical system the ray is deflected by 
a definite angle in one of a number directions depending on the inclination 
of the mirror the ray was reflected by. Change of the direction of the ray 
during multiple passage through the system have therefore the character 
of random walk and may be described mathematically by two-dimensional 
diffusion equation. Finiteness of the aperture of the optical system may 
be approximately taken into account by null boundary conditions at the 
boundary of some finite region of the plane. The region is a disk if the optical 
system has axial symmetry. In this case the problem may be easily solved by 
the method of separation of variables so that and explicit expression in the 
form of an infinite sum may be obtained for the resulting angular distribution. 

Moreover, since each term of the sum exponentially decays with the num- 
ber of passages increasing and the exponents are different, only one of the 
terms (having the minimal exponent) dominates in the case of a large number 
of passages. Thus, asymptotic (corresponding to large number of passages) 
angular distribution of the radiation turns out to be universal and is described 
by the Bessel function of order zero. 
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The reflection of tlie ray by tlie mirrors is considered in tlie framework 
of geometric optics and therefore the method is apphcable only in the case 
of large Fresnel number, a} j XL ^ 1 (here A is the maximal wave length in 
the wave packet, a the radius of the mirror and L the distance between the 
multilobe mirror and the opposite plane mirror closing the optical system). 
This means that the ratio A/a should be much smaller than the ratio a/L. 
We assume too that the thickness of the multilobe mirror h <^ L. 



2 Double passage of a ray through the optical 
system 

If a ray of light which is directed along the unit vector k falls onto the plane 
mirror having the unit normal vector n, then the reflected ray is directed 
along the unit vector (Fig. |l] b) 

ki = k-2(kn)n (1) 




Figure 1: An optical system equipped by a multilobe mirror: a) The general 
scheme; b) the form of the ray reflected by a lobe of the multilobe mirror 
and then by the closing mirror. 

Let no be the unit vector along the axis of the optical system. Assume 
that the vector k characterizing the ray is close to no and the normal n to 
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the mirror is close to —no- Then we can represent these vectors in the form 

k = no + a, n = -no + 7 (2) 

where the vectors a and 7 are orthogonal to no and the terms of the second 
order in a and 7 are omitted. The absolute values of the vectors a and 7 
are equal to the angles of inclination of the ray and the mirror normal to 
the optical axis. The directions of these vectors show in which direction the 
ray and the mirror are inclined. Call therefore the vectors a and 7 'angle 
vectors' or simply 'angles'. 

We shall assume that the inclination angles are small so that only first 
order terms may be conserved in calculations. In this approximation kn = 
— 1 and therefore the reflection law ([l|) takes the form 

ki = -no + a + 27. (3) 

Let now the reflected ray (having the direction ki) be again reflected by 
the mirror closing the optical system from the opposite side. If this mirror 
is orthogonal to the optical axis, i.e. it has normal no, then the direction of 
the ray after this second reflection is characterized by the vector 

k' = ki - 2(kino)no = no + a + 27. (4) 

Finally we see that the ray characterized by the vector angle a converts after 
double passage through the optical system to the ray corresponding to the 
new vector angle a': 

a ^ a' = a + 2j. (5) 

Here the vector angle 7 characterizes inclination of the mirror which reflected 
the ray. This mirror is only one lobe of the multilobe mirror. In each double 
passage the ray is reflected from one of the lobes chosen randomly. 

3 'Diffusion' of the inclination angle and the 
angular distribution of rays. 

Assume that the multilobe mirror contains s mirrors (lobes) characterized by 
the vector angles {71, 72, • • • , %}• After each double passage the ray may be 
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reflected by one of tliese mirrors. If tlie mirror witli number i, i = 1,2, . . . , s 
is ciiosen, tlie direction of tlie ray is clianged as follows: 



d ^ a' — a + 2ji, (6) 

In each double passage the number i is chosen randomly according some 
probability distribution. In the simplest case the probabilities of all lobes 
are equal, Pi = 1/s. We shall consider only this case. 

After double passages an initial inclination angle changes as 

d ^ d' — d + 7) (7) 

where 

^=6 + 6 + ... + 6v 

is a random variable equal to the sum of N random variables ^j. All variables 
have equal distributions. Each of them takes the values {271, 272, . . . , 275} 
with corresponding probabilities. If the probabilities of all values are equal to 
each other, pi — 1/s, then the expectation values, variances and covariations 

— * 

for the components of each of the (vector) random variables are 

= - ^ 27,,, = - ^ 27,,, 

s ■ s ■ 

s ■ s ■ 

* i 

We shall assume for the aim of simplicity that the set of vectors {71} is 
symmetrical so that 

= M4 = 0, Di, = Diy = a\ cov(e,, iy) = 0. 

Then the x- and y- components of the random variables have the same 
characteristics and are independent from each other. For the random variable 
fj (which is the sum of ^i) we have 

M77 = 0, Dr)^ = Dr]y ^ N a"^ cov(77„ rjy) = 0. 
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Change of the angle described as in (|^ is nothing else than Brownian 
walk of a point in a plane. If N is large, then, owing to the central limiting 
theorem of the probability theory, the distribution in various values r = (x, y) 
of the random variable ff may be approximated by the normal distribution: 

Pn{x, y)dxdy = ^^^^^ exp ^ dxdy. (8) 

The probability distribution for the components of the vector angle a' (for 
the given initial angle a) is found readily. If we denote these components by 
the same letters x, y, then this probability has the form 

PN{x,y\a)dxdy = ^-^exp ( 1 dxdy. 

If an arbitrary initial distribution of the angles a is given, then the final 
distribution is 

PN{x,y) = J PN{x,y\d) P{d)da^day. 

If Ncr'^ = tis considered as time, then the distribution function P{t, x, y) = 
PN{x,y) satisfies the diffusion equation: 

dP^lJd^ 

dt 2\dx^ dy^ J ' ^ ^ 

In the preceding consideration is has been assumed that the angle a' may 
become, as a result of the random process, arbitrarily large. In reality it is 
restricted by finiteness of the aperture. This may be taken into account by 
the corresponding boundary conditions for the equation (P). 



4 Angular distribution for finite aperture 

When after a number of double passages the inclination angle becomes suf- 
ficiently large, the ray will be lost owing to the finiteness of aperture. In the 
language of Brownian walk this means that the point disappears and takes 
no part in further random process. If we have initially many Brownian par- 
ticles, part of them is lost because of finite aperture so that the final number 
is less. In terms of the probability distribution, its norm decreases with time. 
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At first glance, this argument leads to null conditions at the boundary 
of some finite region in the plane of Brownian walk. Adding this boundary 
condition to the equation we could find the distribution taking into 
account finiteness of the aperture. In the case of axial symmetry of the 
optical system the permitted region might be a disc of the radius restricting 
the maximum absolute value of the angle depending on the aperture. One 
is attempted to impose null conditions at the angle equal toT = a/L where 
a is the radius of the mirror and L distance between the opposite mirrors of 
the optical system. 

This is however not so simple. Precisely speaking, no null boundary 
condition may be formulated for the angle a' . Yet we shall show that this 
may be made approximately. The null initial condition must be imposed 
at some critical angle depending on F and the angle a characterizing 
deflection at each double passage. 

The precise formulation of the problem is possible only in terms of the 
joint distribution P(a, p) in the angle a of inclination of the ray and position 
p of the point of reflection this ray at the mirror. It must be required that 
this distribution function be zero when p takes values corresponding to the 
edge of the mirror. The problem is that the equation for the function P{a, p) 
is complicated and cannot be dealt efficiently. 

Instead of this correct consideration we shall restrict ourselves by the 
function P{a) and the diffusion equation for it. As for the boundary con- 
ditions, we shall find, up to the order of magnitude, the angle Fc such that 
when this angle of inclination is reached, the ray is with great probability 
lost in some time after this. 

To find the critical angle Fc, we shall argue in the following way. Let the 
inclination angle Fq is reached in the course of Brownian walk. If this angle 
leaves unchanged during next passages through the system, then after each 
double passage the point of the refiection at the mirror is replaced by 2LTq. 
Then after A^'o = F/Fq double passages the point of refiection is replaced by 
the value 2LA'"oFo = 2LF = 2a so that the ray will certainly escape from the 
optical system. 

However after reaching the value Fq the inclination angle will not stay 
constant but change in the course of Brownian walk. The preceding con- 
clusion that the ray will be lost in Nq double passages will be valid if the 
change of the angle is much less than Fq. Making use of typical Brownian 
replacement \/l%a, we have inequality \/I%a <^ Fq as the condition that 
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the ray will be lost. Equivalently this condition may be written as Fq ^ Tc 
where 

r,= (ra2)V3 

Vice versa, if Tq <^ r,, then the random changes of the angle Tq in the course 
of Brownian walk will radically change it before the ray will be lost. In this 
case the fact that the angle reached the value Fq has no special significance. 

We see that the angle Fc is critical in the sense that reaching this angle in 
the course of Brownian walk leads with great probability to loss of the ray. 
This means that the distribution function P{a) must be small for \a\ = Fc. 
Approximately we may take it to be zero for such arguments. This means 
imposing null boundary conditions at the boundary of the disc of radius Tc- 

Thus, the diffusion equation (j^) should be solved with the null boundary 
conditions at the boundary of the disc of radius F^: 

P(t,x,y) = when + = F^. 

This solution is easily obtained with the help of the separation of variables 
in the polar coordinates. It has the form 

P{t,r,ip) = P{t,x,y) 

= XI H 6 * Jn{kmnr) CXmn COS flip 

n=0 m=l 

+ sinmp (10) 

n=l m=l 

where the following notations are used: 

^mn = 2 ^'"^ ~ fJ'mn/^c, ^'mn are the roots of the Bessel function Jn{r). 

The coefficients amn, Pmn are determined by initial conditions as follows: 
P{0,r,ip) = 'Y {amn cos rup + Prnn sin rnp)Jn{kmnr). 

n=0 m=l 

Asymptotically (at large t = Na"^) the first term in the formula (p!0| ) 
dominates because it corresponds to the minimal root of the Bessel functions 
and therefore minimal exponent. Therefore after a large number of double 
passages 

p{t,r,ip) = P{t,x,y) = PN{x,y) ^ aioe"^''^' Ml^ior/T,). (11) 
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It is only numerical factor aio which depends on initial conditions in this 
expression. Therefore, asymptotic (for large A^) angular distribution of the 
radiation issued by the optical system with the multilobe mirror does not 
depend on initial conditions and is described by the Bessel function of zero 
order .[] 



5 Conclusion 

In the present paper we have considered (in the framework of the geometric 
optics) the rays of light in an optical system having an usual plane mirror 
from one side and a multilobe mirror from another side. We evaluated the 
probability that the outcoming ray has a definite inclination angle. It was 
shown that the reflection of a ray from one of many lobes of the multilobe 
mirror inclines it in one of many directions. As a result, the inclination angle 
undergoes Brownian walk. The resulting probability distribution satisfies 
then the diffusion equation, with the number of passages playing the role of 
time. 

If there are many input rays, then the same function which describes the 
probability distribution gives the angular distribution of the outcoming rays. 
It was shown that the distribution alters when the number of passages of 
the ray through the optical system increases. Asymptotically for very large 
number of passages the distribution tends to some universal distribution 
which possesses the same symmetry as the lobes of the multilobe mirror 
(axial symmetry in the considered case). 

If a (modeless) laser is constructed on the basis of such an optical system, 
its radiation will have distribution identical with the asymptotic distribution 
we have found. The angular distribution of the multilobe mirror laser proves 



to be stable under perturbations, therefore the distribution (11) maintains in 
the non-ideal conditions (for example with inhomogeneous active medium). 
This explains why the multilobe mirror laser have smooth (homogeneous) 
angular distribution 0, ^ . 

Several methods could be applied to increase the homogeneity of active 
multilobe mirror lasers. Among them 1) the increasing of the Fresnel number 

^In the very special case when initial conditions correspond precisely zero coefScient 
aio, the asymptotic distribution is determined by the next term of the sum corresponding 
to the larger root of the Bessel function. 
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(i. e. the increasing of the laser aperture and the decreasing of the cavity 
length and wavelength), 2) the increasing of the number of lobes (planes) 
of the multiloba mirror, 3) the application of the active medium with suffi- 
cient broad bandwidth (for example Neodymium glass active elements). The 
above developed method of calculations may be generalized to include these 
modifications increasing the homogeneity of a laser beam. 
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